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CHAPTER I 
INTRODUCTION 

The radar cross section (RCS) of jet intakes has been extens¬ 
ively studied. There are several reasons that it is of interest. 

One is simply that it is a major element in the radar cross section 
of aircraft and must be accurately evaluated to estimate total air¬ 
craft RCS. Another reason for study is that potentially the RCS 
of the aircraft could possibly be reduced if the scattering mech¬ 
anisms are well understood. A third (though certainly not final) 
reason is that many aircraft identification or classification tech¬ 
niques propose to use modulation of the radar return imposed by the 
aircraft engine as a significant identifiable feature. It would 
be a questionable approach to establish such a system on an effect 
which is not well understood, particularly in terms of establishing 
the system's susceptibility to intentional confusion or camouflage. 

This study is somewhat limited to the region for which IXr-X; 
the asymptotic forms developed for higher frequencies may fail in 
this region, and the very low frequency techniques (for which little 
or no energy penetrates the intake) are invalid. 

This study develops two of the practical problems involved 
in calculating the RCS. The coupling coefficients at the mouth of 
the intake are known in princple, but difficult to compute in prac¬ 
tice. The effect of the engine structure is known only generally, 
and has been very loosely approximated in past studies. 

The jet intake can generally be modeled as an open ended wave¬ 
guide with an obstacle (the engine) some distance down the waveguide. 
Figure 1-1 illustrates the geometry. The problem can then, in 
principle, be solved by the generalized scattering matrix technique. 
The scattering matrices of the significant scatterers, if known, 
can be self-consistently manipulated to produce the backscattered 
field. Let 


u 1 represent the incident field, 
u bs = the backscattered field, 

[S,,j = a matrix, representing (in some sense) the direct 
1 J backscatter of the open waveguide 
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Si 2 = radiation characteristics of waveguide modes 

S 2 } = coupling matrix of incident field to waveguide modes 

■ 

S 22 J * reflection of waveguide modes from the open end 

sj = reflection of waveguide modes from the obstacle 
« « 

T 2 b_ = transmission down waveguide 

T b 2 = transmission back from waveguide to aperture. 



Note that off-diagonal terms in each transmission and reflection 
matrix will represent mode conversion from one mode to another. 


From Figure 1-1 we can easily show that: 

bs r i . r 
u = S-j^u S 1 2 r 2 

r l = S 21 u1 + S 22 r 2 
r 2 = T b2 S b T 2b r l 


r 2 = T b2 S b T 2b (S 21 ul + S 22 r 2^ 
[ I_T b2 S b T 2b S 22j r 2 = T b2 S b T 2b S 21 u1 
[ (T b2 S b T 2b^ 1 " S 22^ r 2 = S 21 ul 
r 2 = [ (T b2 S b T 2b ) 1 " S 22] 1 S 21 u 

ubS = { S ll + S 12 t T b2 S b T 2b^ 1 ' S 22] 


( 1 - 1 ) 

d-2) 

(1-3) 

d-4) 

(1-5) 

( 1 - 6 ) 

d-7) 

d-8) 


The matrices S-,p S, 2 » S-i, and S 22 ^ ave been solved for by the Wiener- 
Hopf technique. These results are quite complicated algebraically. 

The heart of the problem, however, is calculation of the Wiener-Hopf 
factorization functions. Considerable effort has been given to cal¬ 
culating these functions, which is presented in Appendix A. This 
study has emphasized circular waveguides, because much prior work 
has been done in the area, and the symmetry of the geometry simplifies 
the problem somewhat. Also, since engine geometries are circular, 
non-circular inlets require an additional model of the mode conver¬ 
sion as the energy travels down a waveguide of varying cross-section. 
Thus our matrices T. ? and T~. are diagonal, having only the relevant 
phase delay for eacfrmode. ^This study does not address further 
development of this effect for non-uniform waveguides. 
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Summary of Circular Waveguide Modes 


: i 

i 

i. 

■i 


Propagation in a circular waveguide is limited to a discrete 
set of modes, which can propagate only for waveguide diameters larger 
than a certain minimum (cutoff) diameter. The notation used in this 
study is consistent with Harrington: TEf 7 refers to a mode with 
electric field transverse to the axis of propagation. In the Russian 
literature (Weinstein) this is referred to as a magnetic mode, since 
it has a z-directed magnetic field, and all other field components 
can be simply derived from it. 

The number of modes that can propagate in a circular waveguide 
goes roughly as the square of the diameter. Table 1-1 lists the first 

TABLE 1-1 



SUMMARY OF LOW-ORDER CIRCULAR WAVEGUIDE 
MODES AND CUTOFF FREQUENCIES 


Mode 

cutoff ^ 

cutoff ka 

TE11 

.5861 

1.8412 

TM01 

.7655 

2.4048 

TE21 

.9722 

3.0542 

TM11 

1.2197 

3.8317 

TE01 

1.2197 

3.8317 


1.3373 

4.2012 

TM21 

1.6347 

5.1356 

TE41 

1.6926 

5.3176 

TE 12 

1.6970 

5.3314 

TM02 

1.7571 

5.5201 

TM31 

2.0309 

6.3802 

TE51 

2.0421 

6.4156 

TE22 

2.1346 

6.7061 

TE02 

2.2331 

7.0156 

TM12 

2.2331 

7.0156 

TE61 

2.3877 

7.5013 

TM41 

2.4154 

7.5883 

TE32 

2.5513 

8.0152 

TM 22 

2.6793 

8.4172 

TE13 

2.7172 

8.5363 

TE71 

2.7304 

8.5778 

TM03 

2.7546 

8.6537 

TM51 

2.7921 

8.7715 

TE42 

2.9547 

9.2824 

TE81 

3.0709 

9.6474 

TM32 

3.1070 

9.7610 

TM61 

3.1628 

9.9361 

TE23 

3.1734 

9.9695 
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28 modes in order of increasing cutoff diameter. The cutoff is cal¬ 
culated by noting the J„(j' )=0 or J'(j' )=0 and then ka=i or 

J nm* $ lnce ka = ^ -^ = it ^ the in terms of ka can be simply 

derived. It seems more intuitive to discuss D/X. Figure 1-3 presents 
electric field line pictures for the first few modes. 

This study employs the time convention; the choice is 
necessary for consistency with the vast majority of Weiner-Hopf 
literature (Einarsson et al). The coordinate system chosen will 
be standard spherical and/or cylindrical coordinates, with the 
origin at the center of the waveguide mouth. See Figure 1-2. The 
polar angle is 0, the azimuthal angle $. The radius will be desig¬ 
nated by r, in spherical coordinates, and p in cylindrical coordin¬ 
ates. The waveguide radius is a, and the waveguide extends from 
z=0 to z=-°° at p=a. 

Since a large part of the literature is concerned with the 
Weiner-Hopf solution to scattering by a semi-infinite circular 
cylinder, an heuristic description will be presented here. Tutorial 
discussions of the Weiner-Hopf technique can be found in [Noble (1958)] 
and [Morse and Feshback (1953)]. 

The Weiner-Hopf technique is based on the fact that the fourier 
transform of a causal function is entire in a half plane. For example, 
the function 

x(t) =|e" at t > 0 

|o t < 0 (1-9) 

has the Fourier transform 


X(o>) = / e 1ut e" at u(t) dt = / e (la) " a)t dt 
-00 o 

= TT^aJ [ e ' a °° - e °] (1-10) 


[Recall that we are using e~ 1ut time dependence]. Hence X(u)) has 
a single pole at to = - ia and is entire (has no poles) for lm(u))>0>-ia. 
In the Weiner-Hopf technique, the fourier transform of the field 
component is taken with respect to the coordinate along one axis 
of the problem (the axis parallel to the semi-infinite object). 

The incident and scattered field are then related by means of the 
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Figure 1-2. Coordinate system for modeling engine inlet. 
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TE 2,1 TE 2,2 TE 2,3 

Figure 1-3. Electric field lines for TE modes 
in circular waveguide. 
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Figure 1-3 (continued). Electric field lines for 
TM modes in circular waveguide. 








physics of the problem. Noble cites three different techniques 
for doing this: field matching. Green's function integral equation, 
and dual integral equations. The resultant equation is manipulated 
so that the left-hand side must be entire in one half-plane and the 
right-hand side in the other half-plane, with a region of overlap. 
These are then equated to a third function, which must be entire 
over the whole a-plane. Because of Liouvilie's Theorem, this new 
function ends up being a constant, equal to zero. This separately 
sets each side of the equation to zero, which leads to the solution 
(by taking the inverse fourier transform). 

For the cylinder scattering problem, the incident field is 
decomposed into cylinder waves by the Bessel function addition 
theorem. The transform is taken with respect to the z-axis: 


e(a,r,<j>) = f E(z,r,<}))e“ 10lZ dz . (1-11) 

-oo 


For each cylinder wave, the incident and scattered fields are related 
by imposing the boundary conditions that the tangential electric 
field must vanish at the cylinder walls, and the surface current 
must vanish in free space. These are used to manipulate the problem 
into a Weiner-Hopf form. This process is described by [Weinstein 
(19C9) and Einarsson et al (1966)]. 

The analysis of interactions at the mouth of an open waveguide 
goes back to [Chu (1940)]. By using the Kirchoff approximation the 
radiation fields for the lowest order propagating modes for a cir¬ 
cular, and a rectangular waveguide are calculated. In terms of the 
GSMT, this would enable us to compute elements of [Si 2]. 

[Levine and Schwinger (1948)] evaluate the acoustic radiation 
and reflection characteristics of a hollow circular pipe. An 
integral equation for the velocity potential is solved via the Weiner- 
Hopf technique. Their study is confined to symmetrical modes incident 
on the open end of the pipe. In a circular acoustic waveguide, 
the symmetrical modes are given by 


*(P*4>) = J o (j 0,m p/a) or J o (j ™ p/a) (1 " 12) 

j_ are the zeroes of J (x) 

om o 

are the zeroes of Jg(x) 


I 


where 




[Pearson (1953)] was the first to apply the Weiner-Hopf technique 
to the electromagnetic problem. He considers a transverse magnetic 
plane wave incident on the open end, and obtains equations for the 
Laplace transforms of the axial and azimuthal currents on the waveguide 
walls. These equations are then solved by the Weiner-Hopf techniques, 
and the currents can be obtained by inverse Laplace transform. The 
fields far from the mouth down the pipe are then evaluated asymp¬ 
totically and found to result from propagating modes, thus giving 
the coupling coefficients. The backscattered fields are not evalu¬ 
ated, nor is the behavior near the mouth studied. 

[Jones (1955)] analyzes the scattering of sound waves by a 
solid semi-infinite cylinder. He considers both hard and soft boundary 
conditions although numerical results are evaluated for only the 
hard (du/dn=0) case. An approximation for high frequency (large 
diameter) is developed, and a low frequency expression is presented. 

A variational expression is developed to establish limits of error 
on the approximations used in evaluating the exact expressions. 

Both the pressure field on the cylinder, and the scattered far field 
are evaluated. The end cap pressure time response due to an incident 
unit step is evaluated by taking the inverse Laplace transform of 
the frequency domain response. 

[Noble (1958)] treats the scalar problem for both radiation of 
the lowest order mode, and coupling of an incident plane wave to 
the lowest order mode in a circular waveguide. His discussion is 
tutorial in nature, drawing somewhat from Jones' work. There is 
a more complete treatment of the general technique used to solve 
Weiner-Hopf problems than is found in most other references. 

[Jones (1964)] applies the Weiner-Hopf technique to radiation 
from a semi-infinite hollow pipe. This study is somewhat tutorial 
in nature, being part of a textbook, and considers only a few of 
the lowest order modes. Jones notes that TE modes radiate more 
efficiently than TM modes. TM modes have higher reflection coef¬ 
ficients at the open end of the pipe. 

[Einarsson et al (1966)] give an exhaustive study of diffraction 
by both the infinite and semi-infinite circular cylinder. The back- 
scatter for a plane electromagnetic wave incident on a solid, semi¬ 
infinite, perfectly conducting rod is given. The scalar (sound) 
scattering from a semi-infinite (both solid and thin walled) tube 
is evaluated. The radiation and reflection for a propagating scalar 
wave incident on an open end are evaluated (reproducing the results 
of Levine and Schuinger and also Weinstein). There is a brief dis¬ 
cussion of finite cylindrical resonators with one end open, one end 
closed (rigid, Dirchlet boundary condition). 
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The general solution for scattering of a plane electromagnetic 
wave from a semi-infinite thin walled* perfectly conducting tube 
occupies about half the report. This includes both the backscatter 
far-field [S]]] and coupling coefficients, for the general plane 
wave (neither TE nor TM). The special case of axial incidence is 
considered. Radiation from a source inside the tube [S 21 ] and re¬ 
flection from an open end CSp^3» are evaluated, largely copied from 
Weinstein. Asymptotic radiation of the far-fields is compared to 
the Kirchoff approximation. 

A substantial part of this study concerns evaulation of the 
Weiner-Hopf factorization functions. A number of forms of integral 
expressions which exactly define them are developed. Power series 
approximations to the low frequency are developed. Unfortunately 
these expressions are still quite complex. 

In addition, only the magnitude of the functions is approxi¬ 
mated; the phase is not. A large effort to develop high frequency 
approximation yields some useful simplifications, but the results 
are expressed in terms of another unknown function, albeit much 
simpler. Numerical data is given for varying values of ka with the a 
parameter fixed. This,is inconvenient since, generally, a=kcos6, 
and we are primarily interested in fixing ka and varying a. Experi¬ 
mental data are presented for finite cylinders. No experimental 
data are presented for a semi-infinite cylinder. 

[Witt and Price (1968)] analyze the problem of a finite tube 
without recourse to Weiner-Hopf techniques. Instead, the direct 
backscatter from the rim, and the coupling coefficients to waveguide 
modes are calculated by what amounts to the Kirchoff approximation. 
The incident field tangential to the aperture is expanded as a sum 
of waveguide modes, and a waveguide admittance for each mode is com¬ 
puted. The reradiation is also calculated via the Kirchoff method, 
using the Stratton-Chu integral. The termination is modeled as an 
impedance, which results in a simple (scalar) reflection coefficient. 
The reflection coefficient is then transformed to its equivalent 
impedance as seen at the mouth of the waveguide. The total waveguide 
admittance for that mode is then calculated, and the scattered field 
is expressed as a sum over the waveguide modes, times the incident 
field projection on that mode, times the generalized admittance. 

For backscatter with vertical polarization (TE) (which corres¬ 
ponds to our $ polarization) they present the formula 


Ey(x' ,0,z') = - M cos . 8 e -jkr 


(1-13) 


vJ^v+Vq) 


V v+ V 


l I E 


n=0 m=l 


nm 


(v)E__(V ) (1+r ) /l + 

nm 0 nm ' V y / 

\ ’nm / 





where 

v„ = kasin0 (incident) 

o o 

v = kasine (scattered) 

E = weighting coefficient for nm-th mode 

nm aa 

Y nm = propagation coefficient for nm-th mode 

r = reflection coefficient seen at aperture for nm-th mode, 
nm 

The work of Weinstein* in many cases predates the works reviewed 
here. However, most of Weinsteins's papers were published in Russian 
journals. His book [Weinstein (1969)] contains virtually all of 
his earlier work, and is readily available. It should simply be 
noted that Weinstein's work did frequently predate work in the West. 

[Weinstein (1969)] has collected all of his earlier work on the 
Wiener-Hopf technique into a single volume. The book is tutorial 
in nature, beginning with the simplest problem, diffraction and 
radiation by a plane parallel plate waveguide. Having developed 
the basic Wiener-Hopf arguments, he proceeds to analyze circular 
waveguides, first considering only symmetrical modes (to eliminate 
azimuthal dependence from the problem). Acoustical problems are 
then solved, followed by the general problem for electromagnetic 
waves scattered and radiated by a semi-infinite circular cylinder, 
including azimuthal dependence. Comparisons with answers obtained 
via the Kirchoff method are frequent, showing those circumstances 
under which the Kirchoff method works and those under which it fails. 
There is a substantial discussion of the relative accuracy of Huygens 
principle, compared with edge diffraction; a substantial point is 
made that edge diffraction yields more reliable answers. 

[Kao (1970)] presents a completely novel approach to scattering 
from cylinders. He determines the currents on finite cylinders by 
using point matching and then calculates radiation patterns from 
the currents. For the semi-infinite cylinder [1970b] he sets up 
two sets of points with slightly different interpoint spacing. By 
various manipulations of these, he determines the magnitude of the 
traveling wave launched on the cylinder by the incident plane wave, 
as well as the current in the vicinity of the aperture. However, 
his analysis is confined to broadside incidence (0=90°). 

[Bowman (1970)] develops ray-optical diffraction expressions 
for the scattering from the end (aperture) of the semi-infinite wave¬ 
guide, and compares these with an asymptotic approximation to the 


*Also transtated as Vajnshteijn, Wainstein, Vainshteir. 







exact Weiner-Hopf solutions. For the asymptotic form of the Wiener- 
Hopf solution, he obtains, for direct backscatter (0=0) 


E BS ^ - x e 


ikr , 1 + e 17T/4 J i m m" 3/2 e i2mka 1 

AFa m= 1 J 


whereas by ray optical methods, he obtains 


(1-14) 


e bs ^ - 


;f F e ikr 


iir/4 00 

5 — l 

Aka m=l 


1 -m e i2mka 

0 m-l m 
Z m 


(1-15) 


) 

For bistatic scattering with axial incidence he obtains from approxi 
mating the Wiener-Hopf solution letting his 0 go to ir-0 and his <j> 
to 2n-<p 


E s x 


, t e ikr ( 2a N^ 2 cos(kasin0-ir/4) (sin<t>) 
+ * ~F~ cos(0/2) 


f, + cos 2 (e/2) ? f 

SWa cose mil 


m -3/2 i2mka 
m e 


From the ray optical approximation he obtains 


(1-16) 


-ikr / \l/2 cos(kasin0-n/4) 

' s -iVfe e ' s,n * 


1 + 


. A/4 

Aka 


cos (0/2) y 
cos0 L 


® -m i2mka 


l e 


m=l 2*^ A 


(1-17) 


- 1/2 

In both cases it is seen that the term of 0(ka ' ) is identical 
up to the first two terms in the summation 


vs 


• e i2ka _ .2 e 4ika _ 1 3 e 16ka 
r~ + 372 ,3/2 


• e i2ka i 2 e 14ka ^giSka 

1 — + ^72 + —f 


(1 -18a) 


(1- 18b) 
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He attributes this difference to the ray-optical approach, for which 
he first considered scattering by plane parallel plates; multiple 
scattering was effected by assuming a cylinder wave was generated 
from each edge after scattering. This result was then specialized 
to having a single point participate in the scattering on each edge 
but not modified to take account of the fact that a cylindrical wave 
is no longer being emanated from each edge. Beyond this comment. 
Bowman does not analyze further, since the exact result from the 
Weiner-Hopf solution instructs us how to modify the ray-optic con¬ 
tribution. 

A comparison of Bowman's results with those of Witt and Price 
is illuminating. If we restrict ourselves to on-axis backscatter 
and disregard terms due to reflection from the termination, the 
formula presented by Witt and Price reduces to 


E 


s 



00 




= -jka 


a 

Jr 






E nm^ 




d-19) 


We observe that, besides the opposite time convention, there is an 
additional factor of jka present in the form presented by Witt and 
Price. They remark that this is recognizable as the physical optics 
approximation for scattering from a flat conducting disk. Since 
this is not in agreement with the high frequency behavior of the 
Weiner-Hopf solution, we conclude that the disk is not a good high 
frequency model for this problem. In fact, it will be seen later 
that the disk does adequately model the on-axis scattering of a 
cylinder terminated by a perfectly conducting plate. 


[Moll and Seecamp (1970)] present a more realistic model of 
the engine geometry. The approach of Witt and Price is used to model 
the scattering, coupling, and radiation at the duct inlet. The study 
is confined to TE modes. The termination, however is modeled by 
two sets of blades, each as shown in Figure 1-4, to simulate the 
first stage of a compressor. The blades are modeled as being planar 
(normal to the z axis). The two sets of blades had different numbers 
of blades and blade widths and assumed varying relative orientations 
(stator to rotor). The scattering at the termination is modeled 
by a similar procedure to that at the inlet. The backscattered field 
is expressed as a sum of modes traveling toward the mouth. The total 
tangential electric fields must match at points where there are 
blades. For each incident mode, integrals were taken over the area 
covered by either set of blades, forcing the fields to vanish. 

The equations thus obtained are used to solve for the scattering 
coefficients for modes generated at the termination. Then, the 
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radiation from each mode is computed by matching the radiated field 
plus the incident field to the internal field at the duct aperture, 
where the internal field is the sum of the modes traveling down the 
waveguide plus those scattered back by the termination. The RCS 
is computed for various orientation of the blade structure and various 
relative orientations of the two bladed structures, giving a range 
of modulation of the RCS caused by the rotor motion. They used 31 
and 37 blades on the respective blade structures. There is a brief 
discussion of a non-planar cap in the duct termination, but the idea 
is not developed. 

[Lee et al (1973)] are primarily concerned with the measurement 
errors made when field strength is measured with a sensor boom. 

The effect of the presence of the boom is analyzed via the Wiener- 
Hopf technique, and the relative distortion thus introduced is cal¬ 
culated. The relevant part of this paper deals with the calculation 
of the Wiener-Hopf factorization functions. Although the general 
factorization expression is developed (as an infinite product of 
factors), the general expression is clearly too complex to be useful. 

A low frequency form is developed for the n=l case for both L+ and 
M+ functions. The n=l case is relevant to low frequencies because 
the TEii mode has the lowest cutoff frequency, hence is also the 
slowest decaying evanescent mode when all modes are cut off. Un¬ 
fortunately the low frequency expression presented do not appear 
to fit very well with data computed in this dissertation by numerical 
integration of the exact integral defining the functions. 

(Lee et al) indicate a constant phase for the low frequency, but 
it appears that the phase varies rather rapidly as the frequency 
increases from O.C. The formulas do work out to the correct "DC" 
form. 


[Mittra et al (1974)] present a detailed report which includes 
reconciliation of the numerical solution with experimental data. 

The complete Wiener-Hopf solution for the semi-infinite cylinder 
is presented. There appear to be some errors in the results presented, 
since the scattering coefficients for direct backscatter (Equation 
(3-86)) could not be reconciled with those presented by [Bowman (1970)]; 
also the two direct polarization solutions for S qq and do not 
reduce to the same answer on axis, nor do they reduce to xhe answer 
given by [Einarsson et el(1966)] for on axis backscatter (Equation 
(5-63)). The solution for large pipes is simplified by an asymp¬ 
totic approximation for the L+ and M+ functions; these are expressed 
as functions of the series 

l m" 3/2 e f"(2ka-n7t + it/2) (1 . 20) 

m=l 


and the same summation, including either even terms or odd terms. 











The Generalized Scattering Matrix Technique (GSMT) is used 
to formulate the problem. Considerable effort is expended expressing 
the individual elements of the matrix in terms of the Wiener-Hoof 
solution. 

Solutions for other geometries are developed by the ray-optical 
method. The geometries considered are an ellipsoid, elliptical plate, 
and semi-infinite elliptical cylinder. The semi-infinite elliptical 
cylinder is analyzed in terms of diffraction by the edge of the 
cylinder. 

These solutions are then combined to estimate the total RCS 
for an aircraft. Numerical results were computed for selected fixed 
frequencies for two aircraft. For calculations involving a termination 
in the jet intake, this was modeled earlier as a perfect conductor, 
or as a dielectric plug of infinite depth, or as a dielectric plug 
of finite depth. 

[Chuang et al (1975)] present an extension of the report by 
[Mittra et al (1973)]. Starting with the exact solution for bistatic 
scattering from a semi-infinite pipe, they develop approximate high 
frequency expressions for the factorization function, based on an 
asymptotic evaluation of the integral defining these functions. 

Both factorization functions can be rather simply expressed in terms 
of a third function, called a modified Lerch function of order 3/2, 

L(x.v) - I »f v e i2m * (1-21) 

m=l 


They then derive via the Mellin transformation a twelve-term series 
representation with complex coefficients, which enables the compu¬ 
tation of this function for v = 3/2 to be mechanized trivially. 

The results thus obtained show excellent agreement with results ob-' 
tained by direct numerical integration of the defining integral, 
down to the cutoff frequency of the lowest mode of that order. Finally 
there is a simplification of the infinite sum which is present in 
the scattering calculation mode by use of the asymptotic form with 
one of the Bessel function addition theorems. 

There appear to be a few errors in this paper. These are summar¬ 
ized in Appendix B. 

[James and Greene (1978)] indicate that both theoretical and 
experimental results show substantial sensitivity to wall thickness. 
They indicate that the exact Wiener-Hopf solution based on infinitely 
thin walls breaks down when the wall thickness is larger than ,0U. 

The results of their summary seems to be that the radiation patterns 
of thick walled pipes are narrower than given by the Wiener-Hopf 
solution. Beyond noting this effect, we will not further discuss 
this problem, since it would require an entirely separate study. 
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CHAPTER II 

DISCUSSION OF WIENER-HOPF SOLUTION 


The exact Wiener-Hopf solution to the diffraction by a semi¬ 
infinite circular waveguide can be found in [Einarsson et al (1966)]. 
The significant formulas from this report are summarized in Appendix 
B. The major numerical problem — that of calculating the factori¬ 
zation functions -- is discussed at length in Appendix A. This chapter 
discusses the physical significance of some of these formulas. We 
reiterate that the e -ltilt time dependence is assumed and suppressed. 

The coordinate system is a substantial stumbling block since 
the incident and scattered field are defined with respect to 6 and $ 
unit vectors which are themselves functions of angle. Further com¬ 
plicating matters, [Einarsson et al (1966)] and [Mittra et al (1974)] 
use different coordinate systems. This study uses the coordinate 
system illustrated in Figure 2-1, which is the same as [Mittra et 
al (1974)], since it forms a self-consistant reference for scattering 
calculations. 

The incident field is assumed to come from the angle (Mu. 
Looking in along the y-axis toward the origin, we see Figure 2-2. 
Looking in along the x-axis toward the origin we see Figure 2-3. 

Hence we obtain the following unit vectors. 


0 = - x cos6 1 - - 2 sin0.j 


(2-la) 


$ = - y . (2-lb) 

For th? scattered field, we are not necessarily constrained to <|> s -it 
(for bistatic scattering). Hence, we obtain general formulas for 
the unit vectors. 

0 = x cos@ s cos 4 > s + y cos0 $ sin<f> s - z sin© s (2-2a) 


$ = - x s in(j> s + y cos$ s 


(2-2b) 


In our coordinate system, the incident and scattered fields are in 
the same coordinate system. In [Einarsson et al (1966)], the direction 
of the z-axis is reversed, but 0,. and the unit vectors for the inci- 
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dent field are defined in such a way that they physically point in 
the same direction as in our system. Thus cose., sine*, E^, e] are 
all the same in Einarsson's and our coordinate systems. (See Figure 


2-2 in Einarsson; E^ occurs for g=0, E] for g = 7 with a being re- 

placed by 0 .). In the scattered coordinate systems, § and $ are 
reversed in'sign (hence also E 0 and E*), and cose has the opposite 
sign, but sine is unchanged. D In addition <j> is replaced by -a 
since the x-axTs is the same in the two systems. 


A. On-axis Results 


The most elementary case of the Wiener-Hopf solution can be 
found by taking the limit as 0-0. In taking this limit we will 
consider both the backscatter and coupling of energy into the wave¬ 
guide. 


Taking the limit of Equation (B-3) (which gives the general 
bistatic scattering from the rim for 0 incident and scattered 
fields) with < j > = tt , e .=0 =6 (in the limit 6 - 0 ) we obtain 


S 


ee 


where 


1im 4* 

6-0 



e n = Neumann epsilon = l;n=0 

= 2 ;n ?<0 

J () is the cylindrical bessel function 
n P 

k is the wavenumber = -y 

a is the pipe radius. 


of the first kind 


L+, M+ are the factorization functions; these are functions 
of (n,a,k); k is always understood and n is usually 
obvious from the context. Hence the only argument that 
is explicitly given is a. L+(k) means L+(n,a=k,k) 


f n = nL+(k)/2kaM+(k) 



1 im 
6-0 



00 



2 (-l)(ka ) 2 

4 [l+ (k)] 2 
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-5-5 - 5 - 5 - (n=l understood) 
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(2-4c) 


since f is identically 0 . 

When combined with Equation (B-2), we obtain 


ikr 


-ia 2 k 


4k 2 a 2 M+(k) 2 -L+(k) 


2 (n=l understood) (2-5) 


This can be seen to be the same as Equation (5-163) of [Einarsson 
(1966)] with r+z; the same procedure applied to (which gives 
the general bistatic scattering from the rim for jj; incident and scat¬ 
tered fields) leads to exactly the same result. 
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(n=l understood) 


( 2 - 6 d) 


We see immediately from Equations (B-5) and (B- 6 ) (which give 
the cross-polarized backscatter) that sin mr =0 so no cross-polarized 
field is generated. 


The radar cross-section is defined by 
. 2 

o = lim 4irr 

r-00 



(2-7) 


Applying Equation (2-5), we obtain 
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(2-3a) 


(2-8b) 

(2-9) 


The results of Equation (2-9) are evaluated and plotted in Figure 
2-4. It is of interest to note that the normalized cross-section 
peaks at D/X=.55, with a/na^=10.4 dB, and then decreases very rapidly. 
The lowest order propagating mode (TE11) is enabled at D/X=.586; 
hence this peak occurs just below cutoff for this mode, and as the 
mode is able to transport energy, the cross section rapidly decreases. 
The first three peaks and corresponding modes are summarized in Table 


2 - 1 . 


TABLE 2-1 

SUMMARY OF ON-AXIS RCS PEAKS AND MODE ACTIVITY 


Location of 
peak (D/x) 

Height of 
peak (dB) 

Cutoff for 
mode (D/x) 

Mode 

.55 

10.4 

.586 

TE11 

1.68 

3.5 

1.697 

TE12 

2.69 

2.5 

2.717 

TE13 


The rapid variation of cross section in these regions is related 
to the strong coupling of the normally incident plane wave to the 
TEi m modes, as discussed below. 

[Bowman (1970)] and [Chuang et al (1975)] use their asymptotic 
forms for L+ and M+ to obtain a simple form for the on-axis cross 
section. 
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Equation (2-11) is plotted in Figure 2-5 and yields good results 
for D/X>1. 

Taking the inverse Fourier transform of Equation (2-5) with 
the spatial dependence suppressed, we obtain the time domain response, 
shown in Figure 2-6. This represents the backscattered field as a 
time function resulting from a normally incident plane wave, impulsive 
in time. This result was generated via the discrete Fourier trans¬ 
form of the backscattered formula for 0<D/A<3.2. Naturally the time 
domain response is dominated by the peak at D/A=.55, which results 
in damped oscillations. The time domain response can be related 
(perhaps dubiously at low frequencies) to a ray optic model of scat¬ 
tering at the mouth. The optic model was developed by [Bowman (1970)]. 
For example, if we take the simplest asymptotic approximation for 
the L+ and M+ functions, given by Equations (A-27), (A-29), (A-31), 
(A-58) and (A-59), and substitute them into Equation (2-5), suppres¬ 
sing the spatial decay and propagation, 
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retaining only terms of 0(1) and 0(ka ), with a unit incident 
field 
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Taking the first few terms of the summation explicitly 
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Figure 2-6. Inverse Fourier trans 
backscatter with weight' 
Gibbs-type ri 
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Note that e c results in a time delay of . We employ the Fourier 
transform pair 

1_ ^ Mill 

✓w /2rrt 

and its Hilbert transform 
J_ < > -u(-t) 

✓to /-2ut 

Taking the transform of the frequency domain expression 
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leads to 
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This is plotted in Figure 2-7. 

We note in Figure 2-6 that the sharp peaks at t=m.|^ become 

almost negligible after m=3, and the response shown in the discrete 
Fourier transform results seems to be dominated by the ringing as¬ 
sociated with the lowest frequency resonance. Since this resonance 
is not well predicted by the asymptotic forms (compare Figures 2-4 
and 2-5) it is not surprising that the long-time or steady-state 
response is not well predicted. This resonance appears to be related 
to the currents excited at the mouth, and to some extent, associated 
with an exterior natural resonance of the cylinder, since all the 
waveguide modes are evanescent in this region. 
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: igurf> 2-7. Inverse Fourier transform of asymptotic on-axis backscatter. 











For purposes of comparing the relative importance of these 
coupling coefficients, the power flow associated with each of the 
modes can be computed. For TE modes, with an incident plane wave 
of unit amplitude based on [Collin (1960), p. 179] the power flow 
is computed as 
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We assumed a unit incident field, which has an incident power 

1 Fa i 2 1 Fa 

density of E = so the normalized power flow becomes 

(normalized to the average power flow and the area of the waveguide 
mouth) 
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Similarly, for TM modes, we obtain 



koinm ("T~—) 
nm J™, 
nm 


nm' j ' 
nm 


// |E Z | 2 


(2-22a) 


2tt a 

I I 

O 0 




cos^n 


*l A nml pdpd4> 


(2-22b) 



(2-22c) 


(2-22d) 




( ka >Km a) 

ST 



(2-23) 


The power flow coefficients are plotted in Figure 2-8 for the 
first three TE-j m and TM^ modes, for axially incident plane waves. 

The reason for the behavior of the on-axis cross section becomes 
somewhat clearer, since we see that the TE-, modes all couple much 
more strongly to the axially incident plane m wave than the TM, modes. 
Hence, the region in which they are enabled exhibits a much more 
dramatic variation due to the size of the change of power absorbed. 
[Weinstein (1969), p. 151] discusses this behavior and proves that 
asymptotically the optical cross section equals the sum of the ab¬ 
sorption cross sections for the TE, modes. Figure 2-8 can be com¬ 
pared with Weinstein's Figure 53. 
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Figure 2-8. Coupling of on-axis incident 

modes - normalized p 







B. 


Off-axis Behavior 


The radiation patterns of the various waveguide modes have 
been produced in the literature; all that is produced here are a 
few curves for comparison. Figures 2-9 and 2-10 with ka=4 can be 
compared with [Weinstein (1969)] Figures 46 and 48. Figure 2-11 
and 2-12 with ka=12.77 can be compared with [Narasimhan (1979)]. 

Of more practical importance is the relationship of these 
radiation patterns to radar cross-section. As is discussed in Section 
II-C, the radiation pattern of a single waveguide mode in a certain 
direction can be seen, by reciprocity, to be directly proportional 
to the coupling of an incident plane wave from that direction to 
the appropriate waveguide mode. If, for the moment, we assume that 
the cylinder is terminated with a perfectly conducting flat wall, 
each propagating mode will be reflected back to the waveguide mouth 
unattenuated. For a monostatic radar system, this will mean that 
at some angle, which happens to couple well to a particular waveguide 
mode, that waveguide mode will radiate equally well in that direction. 
Therefore, it becomes of considerable importance to know the directions 
and relative strengths of the coupling coefficients of the various 
waveguide modes. Obviously, for 9=0, we can see from Figure 2-8 
that the TE-j-j mode dominates all others. From Table 1-1 we see that, 
for increasing waveguide diameters, the TMg^. TE^ 7 , TM^, TEgi are 
the next modes to propagate. In a manner analogous to Figure 2-8, 
the relative importance to radar cross-section for these modes is 
shown in Figure 2-13. However, Figure 2-13 differs in that the 
direction is varied for each mode so that the incident field is 
assumed to come from the optimum direction. This direction is plotted 
in Figure 2-14. For example, for D/x=1.5, the TEq 7 has a maximum 
for Ex at 30°, the TMq] and TE 21 modes also have maxima at 30°, but 
for Eq. However, the relative magnitudes are 0.7, -6.8, and 1.0 
dB, respectively. Figure 2-13 might be termed the cavity cross- 
section, since it is the contribution to the RCS which coupling to . 
and radiation from a single mode would produce, assuming that only 
that mode contributed to backscatter. Naturally, this is not true, 
as other modes will contribute to some extent, as well as the direct 
backscatter from the rim. We see that for on-axis scattering, the 
"cavity cross-section" produced by the TEmode is remarkably similar 
to that of the disk. 

For an incident plane wave, we have from Equations (B-9) and 
(B-10) (which give the coupling of incident plane waves to the axial 
waveguide field) 
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2-9. Radiation 
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Figure 2-11. Radiation pattern for TEmode with ka=12.77 (D/X=4.065 
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Figure 2-14. Angle for beam maximum for six lowest-order 

waveguide modes. 
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Fran Equation (B-17), (which give the radiation patterns for the 
waveguide modes) we obtain scattered fields of the form 
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(2-25) 


Consider first the topline of Equation (B-9). 
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This gives us E” m = A^ m () in Equation (B-15) except for the 

direction of propagation. If we assume perfect reflection, and 
ignore the sign change, (since we are interested only in magnitude) 
this leads to 
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The RCS becomes 
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Since we chose 4>..=m, we likewise set 4 > s =tt and vary 6=0 1 -=e s to its 

optimum value. Note that for computations involving A^ m and B 0 m , 

there is a sin n<j> dependence in the modal fields. To adjust this 
to conform with Equations (B-15) and (B-16), it is necessary to rotate 
the coordinate system. This results in radiation patterns evaluated 

at 4> $ * Hence, we obtain, for the remaining three cases: 
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The backscatter directly from the rim is discussed by [Chaung 
et al (1975)]. They used a pulse radar and extracted the RCS of the 
rim based on the first return pulse. The agreement between theory 
and experiment is better in some cases than others, but fails to match 
the detailed pattern. The normalized RCS can be derived from Equations 
(B-l) through (B-4) as follows. 
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Comparison of Figures 2-15 and 2-16 with Figures 2 through 

5 of [Chuang et al (1975)] reveals that despite use of the complex 
scattering form, and more accurate computation of the factorization 
functions, some substantial discrepancies still exist between 
theory and experiment. Possihly the exDlanation is that 
any real cylinder must, of necessity, have a finite wall thickness. 
[James and Greene (1978)] showed that this leads to substantially dif¬ 
ferent results than obtained with the infinitely thin "knife-edge". 

Of course, the accuracy of the measured results is a possible source 


of discrepancy. 
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Figure 2-15. Rim backscatter cross-section for semi-infinite cylinder 

with ka=7.261 (D/X=2.311). 








Figure 2-16. Ri m backscatter 




irocity Considerations 


It is evident from physical considerations that there must 
be some relationship between the coupling of an incident plane wave 
to a waveguide mode, and the waveguide radiation pattern of that 
mode. This will be shown by the reciprocity theorem. 

From [Harrington (1969), p. 116-117], if the current J a produces 
fields (E a ,H a ), and current J b produces fields (E b ,H b ), then 


- #(E a x H b - E b x H a )*ds = J/J(E a -J b - ?-J a )dv (2-36) 
s v 

where the surface and volume are of finite extent. Generally speaking, 
reciprocity is applied to sources and matter of finite extent. In 
this case, however, the matter is of semi-infinite extent. Usually, 
it is shown that the surface integral vanishes far from all matter 
and sources. In this case, we will assume that it vanishes external 
to the waveguide, and thus need only evaluate the 'power flow' 

(E a x H^) inside the waveguide. 

First, we apply reciprocity to determine what source will produce 
a plane wave incident on the waveguide. We assume that there exists 
a current djpole, impulsive in space, of either 9 or $ orientation, 
located by R p = (0=0 O , <|>=tt, r=R 0 ), where R 0 » a, R » X, as shown 
in Figure 2-17. This dipole produces a far-field at the origin of 


(2-37 a) 
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for § or $ source unit vectors, respectively. Hence, to produce 
a 'plane wave' of unit amplitude in the vicinity of the origin, it 
is necessary to set 











incident "plane-wave". 
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Let us first consider TE modes. These have fields given bj 
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where e means the fields propagate toward the waveguide mouth 

-ia' 2 

e means the fields propagate toward - °°. 

Let the b source be located within the guide at z = -L, as shown 
in Figure 2-18. 
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This current source will generate the appropriate TE mode of unit 
amplitude propagating toward the waveguide mouth. It will generate 
the same mode, of equal amplitude, propagating toward the infinite 
recesses of the waveguide. Note that since the leading sign on H p 
and changes with direction of propagation, the discontinuity in 
tangential H exactly matches the surface current. However, E E. 
and H z are all properly continuous. p ^ 

Next, according to Equation (B-10), an incident plane wave 
with ( ) polarization produces a field inside the guide. 
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According to Equation (B-17), the TE mode propagating toward the 
waveguide mouth produces a far-fieldpattern 










Plugging the currents and fields into the reciprocity equation, 
taking the surface and volume indicated in Figure 2-18, we obtain 
the following result. For § incidence, H? has sin n<j> dependence, 
so we take the bottom line of Equations (2-41) and (2-42). 
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We invoke the orthogonality of modes to assert that the sum¬ 
mation reduces to a single term, as this integral vanishes except 
when the mode and current have the same radial wavenumbers and azi¬ 
muthal dependence. When this is satisfied. 
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Since is generated by a field with a sin n<j> dependence, 
rotate the radiated field generated by C q m by tt/ 2. Hence, 
ate with a cos n<J> dependence. 
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This surface integral within the waveguide becomes 
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Then, since each of E ft , Ea, H p , and are based on H , the fact 

lples forces^the integral to z vanish identi- 
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which can be seen to be satisfied by some simple algebraic manipu¬ 
lations. 


For a $ incident plane wave, we note that H? has a cos n<j> depen¬ 
dence, so we take the top lines of Equations (2-41) and (2-42) for 
both the current and waveguide fields. This results in exactly the 


same integrals for both E a -and E^-0 a , with so that 
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which is satisfied identically. 
For TM modes, we begin with 
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This will generate the appropriate field traveling toward the wave¬ 
guide mouth. However, to generate the discontinuity in tangential 
H, the leading sign on all terms must be reversed for the mode 
traveling toward the infinite recesses of the tube. However, as 
before, this is of little importance since it contributes nothing 
to the reciprocity integral. 

The incident plane wave generates fields 
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For the 0 incident field we obtain 
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Finally, for $ incidence, the field generated has sin n<j> depen¬ 
dence, forcing us to the bottom line of the current and field distri¬ 
butions (Equations (2-52) and (2-53)). The integrals all end up the 
same, resulting in 
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Therefore, we have demonstrated the relationship between coupling 
of incident plane waves to waveguide modes and the radiation patterns 
of waveguide modes. 
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CHAPTER III 

SCATTERING BY A REPRESENTATIVE ENGINE-LIKE 
OBSTACLE IN A WAVEGUIDE 


This chapter is concerned with scattering by a simplified model 
representing an engine geometry in the intake duct. The model consists 
of an infinite circular waveguide housing an axially symmetric cone 
centered on a flat plate, both perfectly conducting, as shown in 
Figure 3-1. This model is based to some extent on observation of 
the Pratt-Whitney J-57 turbofan jet engine, on display at the Air 
Force Museum, Wright-Patterson AFB, Ohio. In this study, it is 
assumed that the TE^ mode is incident; similar procedures could 
be followed for any other incident mode. The coordinate system is 
centered at the base of the cone, as shown in Figure 3-2. 

The calculation of fields scattered by the cone in situ would 
be known completely if one knew the exact currents on the obstacle. 

The scattered fields can be obtained straightforwardly from the dyadic 
electric Green's function and the dyadic magnetic Green's function, 
based on the integrals given by [Tai (1973), p. 9]. 


E(R) = iu)y o /// § el (R,R')-W) dv' (3-1) 

™ = ;;; 5 m 2 (i ?,w).w) dv . (3-2). 


Although the source singular term for the dyadic electric Green's 
function is still being discussed in the literature by[Yaghjian (1980)] 
and others, the Green's functions at points not near the source can 
be obtained unambiguously from residue series expansions. These 
expansions consist of a double summation over n(azimuthal index) 
and m(radial index) of terms representing TE and TM fields, both 
propagating and evanescent. Hence the coupling to a given mode can 
be determined by a single integral, allowing one to obtain either 
E z or H z , based on known currents. Mathematically this can be ex¬ 
pressed as follows: Omitting the source singular term, one can write 
(using Tai's notation) 
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Figure 3-1. Basic geometry for engine scattering model 



Figure 3-2. Coordinates system and dimensions for 
scattering from an axial conducting cone. 


<c 










G e (*.R') = ii J J V 

n=u m=i 


.2. \ N ®nX^- k X^ N ®n (+k X> + 
l k,I, o o 


'X*X 


♦ i ' -Be. (+k )We* (?k )t 

Tkl o " 11 ' “ o nv " 

* V V j 


for z £z' 


(3-3) 


G^OU* 1 ) = 7 x § el (R,*') 


= I I e n ‘ 

41T n=0 m=l n 


,2. T R ®nX ( t k X^ N ®nX ( +k x) 

X k x I x o o 


Ne„ (+k )Ne' (+k ) 
v\i v o^- » n"^ » 


for z £ z' 


(3-4) 


where 


“ * J™ /a 


X ■ V 3 


k„ • ,/kV * a 


I 

nm 


I„ = ^a 2 (l-n 2 /j' 2 )J n (j ni i) 2 
p nm n' nm 


k x = / k 2 -X 2 = a, 


I, = 

X n nm 


nm 

2 


S' k 


"x * k 


e n = 


n=0 

n=l,2,3. 


7T(h) represents TM electric fields or TE magnetic fields; Jf contains 
an axial component 

fl(h) represents TE electric fields or TM magnetic fields so that 
is purely transverse. 
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For example, for a single TE mode, we can obtain 
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and similarly calculate E using §,. The elements of 5 , and G ? 
are written out in Appendix C. 


Because of the difficulty of and restrictions on previous 
solutions for current, a new technique of solving for currents was 
explored. A few cases have been previously solved with great effort. 

For example, [Tesche (1972)] solved the problem of a skewed wire between 
parallel plates by using images. [Wang (1978)] used the dyadic elec¬ 
tric Green's function to solve for the currents in an arbitrarily 
shaped dielectric body inside a rectangular waveguide. [Harrington 
(1961), pp. 402-406] gives a variational technique to find stationary 
formulas for scattering using approximate current distributions. 

He applies this to a post in a parallel-plate guide. Unfortunately, 
this method requires computing the self-reaction of the assumed 
current distribution, which in most cases is equally as difficult 
as solving the problem exactly, since it requires computing the 
electric field generated by the assumed current in the source region. 
Hence, no computational advantage is obtained over solving directly 
for the actual current by the method of moments, or some similar 
technique. 


In principle, the problem can be solved by the method of moments 
by assuming the current distribution to be a collection of a series- 
of pulses with unknown weighting coefficients, calculating the elec¬ 
tric fields generated by these pulses, and adusting the coefficients 
such that the tangential electric field vanishes as the surface of 
the obstacle. The problem with this approach, as pointed out by 
is that for coplanar field and source, the dyadic 
s function converges extremely slowly, if at all, 
in addition to the problems involved with the source singular term. 
Convergence of the residue series is enforced by the axial propagation 
ict |z-z' | -a n _|z-z'| 

factor e , which becomes e , for large enough m. 

For coplanar source and field, this factor becomes unity. In Wang's 
case the problem was solved by summing the resultant series (without 
the convergence factor) in closed form. In the case of the circular 
waveguide, we have a Fourier-Bessel series whose summation in closed 
form is not readily apparent. For example, the z component of G j 
is given by 


[Wang (1978)], 
electric Green' 
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Taking asymptotic forms for j and J n ( ), assuming n=l, for coplanar 
source and field point, we obtain 
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It seems that even if p and p 1 are different, this series does not 
converge. The series does converge for real objects after inte¬ 
grating over volume currents of finite extent, since integration 
over p’ results in dividing the asymptotic terms by (m+^JiT, and 
integration over z multiplies the element by the differential dz, 
forcing the coplanar elements to make an infinitesimal contribution. 
The reminder of the integral is well behaved, due to the convergence 
factor for noncoplanar points. 

This extremely cumbersome process seemingly cannot be avoided. 
However, poor convergence in the source region suggests a similar 
problem in free-space scattering, and the superior convergence of 
the magnetic field integral equation (MFIE) over the electric field 
integral equation (EFIE) for most obstacle scattering problems^ 

We therefore proceed, analogous to the MFIE to force H s = n x H xo 

just outside the surface of the conducting body and n x H tota ^=0 
just inside the surface. As with the MFIE, taking the mean at the 
surface results in dividing the current by two. Assume the current 
consists of a collection of pulses of arbitrary weighting 


■’s ’ | V (r - r k )5 k < 3 - 9 > 

wheSVe ti. is a unit vector tangent at the surface at r. , and p(r-r. ) 
is localized near r^. Then K 
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(3-10) 


H S (*) = /// G ni2 (R,R , )-3 s (R 1 ) dv 1 


= I a k // G m2 (R, RI )-u k p(R , - r k ) ds' 
Enforcing the boundary condition 


= n x H total = n x H 1 + n x H s 


n x H 1 =*sJ $ - n x H s 


n x H (If) * | a kP u k' n x M G m2^ R,R ^*\ p ^ R, ' r k^ ds ’ 


(3-11) 


We can then enforce this condition as many times as necessary to 
obtain the required number of equations to solve for the a k 's. The 
advantage of this approach is found only by careful examination of 
the convergence properties of G^ and G m2 . The least convergent 

term in G el goes asymptotically as m sin mx, but the least convergent 
term in G m2 goes as 1 sin mx. After integrating over the surface, 

we obtain S1n m mx or —— m mx , which are convergent. Hence, we can 

work with surface, rather than volume currents, and construct a 
relatively simple code. Additionally, there appears to be no dispute 
concerning the behavior of G 2 in the source region. 


With the geometry shown in Figure 3-2, we define 
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(3-12) 


Assume the current consists of the sum of pulses 
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* s = Z a k^ P(p“Pk^ + k Z b k $ P(P _ P|<) 


(3-13) 


where p(p) is a triangular pulse function. This current distribution 
is shown in Figure 3-3. Then using Equation (3-11), we obtain 


n x H 1 = y a. t%>(P-P.) - n x JJ G m ,(R,R')-t p(E'-£j ds’ 


+ I b k 4>*sP(P-P k ) - n x // 5 m2 (R,R')-$ p(R'-r k ) ds'j 
k * 


(3-14) 

Next, to employ Galerkin's method, we generate 2N equations by multi¬ 
plying through by each basis function and integrating over the surface 
S, to obtain 


u^n x // H i (R)p(R-r ) ds = 


N r 

^ a k ^//P^" r k^^“ r q^ ds ^i *^”^i*^ X J'J*J"/P^“ r q^m2^’^ 1 ^ * 


■ 

tp^'-r^dsds' 


+ k I ] b k^//P (R “ r k^ p ( R " r q^ ds ^i'^“^i‘" x /J’/. ! ’ n ^' r q^m2^ R ’ R ^' 

•$p(R' -"r k )dsdsj 


(3-15)' 


q = 1,2,...N 

1 ■ 1,2 


A A A 


where u-| = $ u^ = t. 

Thus 2N linear equations for the 2N unknown coefficients a. and t> k 
can be obtained, which can be straightforwardly solved by linear 
algebra. 

In order to further speed computation, elimination of unknown 
currents in the backplane can be obtained by using the method of 
images. Instead of the original problem, we introduce an image such 















that the image of the cone has image currents on it, and the image 
incident field is assumed to be present, as shown in Figure 3-4. 

This automatically forces nxE^°^ a ^=0 and n-H^°^ a ^=0 in the backplane, 
hence reducing the amount of coplanar integration necessary. With 
this technique, we obtain a solution for the currents on the cone, 
then integrate these with G m 2 to obtain the H-field in the backplane, 
hence the currents in the backplane. The current distribution thus 
obtained is shown in Figure 3-5. It would be extremely difficult 
to measure this experimentally, and constraints of time and money 
prohibit this verification. 

The fields in the region of the cone-base termination can be 
approximated in the spirit of the WKB approximation, if the cone 
is slender enough, and the cone diameter varies slowly enough with 
axial position. At each axial position on the cone, it is assumed 
that appropriate TE fields exist for an infinite coaxial ljne of 
the same inner diameter. Since E z =0, this guarantees n x E = 0 on 
the cone and outer waveguide walls. The radial wavenumber is then 
computed by solving the characteristic equation 


Jp(ub) Y;(ya) 


Y^yb) j;(ua) = 0 


(3-16) 


for the smallest positive y, with b the inner diameter and a the 
outer diameter. The values of y are plotted in Figure 3-6. Next, 
the scalar wave equation for H, is approximated by assuming that 
y(z) varies slowly enough so that we can neglect the coupling of 
y(z) through the radial function. Suppressing the azimuthal depen¬ 
dence, let H = R(p,z) Z(z). The scalar wave equation 


11_ 

P 3p 


9H z 

3p. 


„2 ? 

Vz + ~T + k H z 
p 8z 


is approximated by setting 


= 0 


(3-17) 


R(p,z) = 0 


(3-18) 


so that the wave equation reduces to 

+ 1 3R + ( 2 n* 

2 7 TZ + (y 

SO 2 0 90 0 J 


Z(zH(k 2 -u 2 )R(e)ZUHR(p) - 0. 


3z 


(3-19) 
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Since it was assumed that R(p) was an appropriate combination of 
Bessel functions, this reduces to 


4 ♦ (kViz - o 

3z^ 


(3-20) 


Since p(z) is a known function, based on the local radius, this 
equation can be integrated numerically to obtain Z(z). The results 
of this numerical integration are given in Figure 3-7. Finally, 

p 

to normalize the radial function properly, J|R(p,z)| da must be 


constant locally at all cross-sectional planes. Subsequently, H p , 
H^, and H can be determined at the inner conductor, and hence tne 
surface currents J t and J , which are shown in Figure 3-8. 


There are several limitations^ this approach. Since H p = 0 
on the inner conductor but H f 0, n • H ^ 0. Maxwell's equations 
are only satisfied approximately, not exactly, since the coupling 
of p(z) through the radial functions was neglected. The behavior 
of the fields near the tip of the cone was approximated very loosely 
since it was assumed that the field consisted only of incident and 
scattered TE^ modes immediately beyond the tip. 
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Figure 3-7. Solution for axial variation of fields for TE,, mode 
in circular waveguide with cone. L/a=2.0, b/a=0.5, ka=z.O. 
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gure 3-8. Currents induced on cone in circular waveguide with 
incident TE-|, mode, l/a=2.0, b/a=0.5, ka=2.0, solution by 
coaxial approximation with slowly varying 
center conductor. 











CHAPTER IV 

CONCLUSIONS AND DISCUSSION 


The problem of predicting the radar cross section of jet intake 
cavities in a spectral region where the cavity aperture is of resonant 
dimension has been approached via the exact Wiener-Hopf solution 
for a semi-infinite waveguide. Exact (numerical integration) RCS 
computations have been presented for the semi-infinite circular wave¬ 
guide for electrical diameters from zero to three wavelengths. These 
computations smoothly join high frequency asymptotic results and 
also serve to define the region of validity for the simpler asymp¬ 
totic models. Therefore one of the matrices in a generalized scat¬ 
tering matrix development of the scattering by a finite loaded circular 
waveguide has been completed. In the process, the nontrivial relation¬ 
ships between two earlier studies of the semi-infinite cylinder have 
been developed and certain errors in a publication from one of these 
studies have been corrected. 

The coupling coefficients relating an incident plane wave to 
internal waveguide modes have been extracted from earlier studies, 
recast in a more convenient form, and evaluated in the resonance 
region. The predominance of the TE modes over the TM modes in power 
absorption for axial incidence has been demonstrated. The relative 
importance of the five lowest order modes for non-axial incidence 
has been demonstrated graphically. The relationship between the 
coupling coefficients and the radiation patterns of the waveguide 
modes has been established via the Lorentz reciprocity theorem. 

The relationship thus demonstrated is that the radiation pattern 
and the coupling coefficient for a given mode at a given frequency 
and polarization are proportional. This is of particular importance 
to RCS calculation since a given waveguide mode will couple to and 
radiate efficiently in the same direction and with the same polari¬ 
zation. 

Modelling of the leading surface of a jet engine as a cone 
on a flat plate inside a circular waveguide, we have computed the 
currents on this obstacle, by using the method of moments and the 
dyadic magnetic Green's function appropriate to the interior of a 
circular waveguide. The reflection coefficient for an incident TE-^ 
mode has been calculated based on these currents at selected fixed 
frequencies. 

It therefore remains to apply these results to Equation 
(1-8) to solve the radar scattering problem in a self-consistent 








fashion. This is done for some illustrative cases. First, fixing 
a perfectly conducting flat plate at ten radii down the waveguide 
from the mouth as shown in Figure 4-1, the RCS has been computed 
while varying the diameter from zero to 1.5 wavelengths. The normal¬ 
ized RCS is shown in Figure 4-2. In Figure 4-3, the normalized RCS 
is shown for the case when the waveguide diameter is one wavelength 
and the flat plate varies from two to ten radii down the waveguide. 

Piscussion 

The calculation of scattering from the inlet mouth and coupling 
to internal modes has been performed modelling the inlet as a circular 
duct with infinitely thin walls and a knife-edge rim. In practice 
inlet rims have a finite wedge angle, and are very rarely circular. 

The major inadequacy encountered in calculating rim scattering for 
the knife edge by ray optics proved to be the underestimate of the 
peak just above cutoff of the lowest order propagating mode (comparing 
Figures 2-4 and 2-5). Unfortunately, for arbitrary geometries, there 
seems to be no simple way to patch up ray-optics to solve this problem 
Since the frequency range this occurs in is low, it may be possible 
to evaluate this region by moment method techniques. Although this 
was not done in this study, we can see that the ray optics approxi¬ 
mation does work well provided we are sufficiently above cutoff of 
the lowest order mode. For circular geometries, this condition is 
that the diameter is greater than one wavelength. 

Because of the reciprocity relations demonstrated in Section 
II-C, it is clear that discussion of radiation patterns and coupling 
coefficients are redundant. The calculation of radiation patterns 
has been discussed by [Weinstein (1969), pp. 139-150]. Weinstein's 
comments on the use of Huygen's principle for calculation of radiation 
patterns indicates that in the transition region (near cutoff) 

Huygen's principle performs poorly, and worse for TM modes than TE 
modes. As Figures 2-13 and 2-14 indicate, the relative importance 
of TM modes to backscattering from cavities is very far below that 
of TE modes, and limited to angle far from the forward direction. 

For all practical purposes, forward scattering in the transition 
region can be calculated using the TE-j-j mode, adding the TE 21 and 
TEqi modes as requirements may dictate. The reason for this can 
be grasped physically by considering Figure 1-3. For modes higher 
than m=l and for all TM modes, the oscillation in sign across the 
waveguide cross-section forces the average interaction to a very 
small value, regardless of incidence angle. For m=l modes, the vari¬ 
ation due to cos n<j> dependence can be seen to lead to an optimum 
angle, from which the aperture fields appear to be nearly of the 
same sign. 

For non-circular geometries, the primary conclusion we might 
draw is that the appropriate low-order TE modes will dominate scat¬ 
tering in the forward direction. Since Wiener-Hopf solutions are 























Figure 4-3. Cross section of semi-infinite circular waveguide 
for axial incidence, D/X=1.0, with flat plate at two to 
ten radii down guide. 









not available for these arbitrary geometries, it will probably be nec¬ 
essary to construct moment method solutions for near cutoff frequen¬ 
cies. Higher frequencies can probably be adequately approximated by 
ray-optic solutions as done by fPathak and Huang (1980)], for example. 

The impulse response waveform given by Figure 2-6 indicates 
that rim backscatter is really only significant for times of t<5a/c. 
Furthermore, this short time response can be reasonably approximated 
ray-optically, confirming our earlier statement that high-frequency 
scattering can be adequately modelled ray-optically. The low- 
frequency scattering, which is not ray-optic, corresponds to the 
long-time response, which is typically of little interest. 

Having enumerated the difficulties in calculating scattering 
coefficients for object inside waveguides, it becomes clear that 
much work remains to be done. The use of the magnetic field dyadic 
Green's function and appropriate boundary conditions provide some 
improvement over use of the dyadic electric Green's function, but 
it is still very expensive computationally. The possibility of 
azimuthal asymmetry (for example, a blade structure) has not yet 
been addressed. One might model the blade structure axisymmetrically 
by using boundary conditions in the backplane E p = 0 and = H p = 0 
and at the outer wall J p (P=a) = - H ( j ) (p=a) = 0. ^This allows radial 
currents and azimuthal E-fields to exist, but not radial E-fields 
or azimuthal currents. It further forces the axial current to vanish 
at the tip of the blades. For a structure with 28 blades, for example, 
it would require an extremely large diameter waveguide before signifi¬ 
cant azimuthal currents could flow on the blades. One problem in 
implementing this study would be that the image procedure used in 
Chapter III to eliminate coplanar integrations could not be used, 
and it would be necessary to integrate coplanar source and field 
points. Another approach would be the use of the free-space Green's 
function, to analyze the structure. In this case it would be 
necessary to also set up currents and enforce boundary conditions 
at the waveguide walls, thus vastly increasing the number of current 
elements required. Another possible approach is to segment the 
obstacle structure and use waveguide modes appropriate to coaxial 
waveguides near the obstacle. However, the convergence properties 
of this approach are unknown. 
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APPENDIX A 

EVALUATION OF WIENER-HOPF FACTORIZATION FUNCTIONS 


The factorization functions L + (oi), and M + (ot), are defined as 
being functions which are analytic in the upper half of the complex 
a-plane, which also satisfy the equations 

L + (ot)L + (-ot) = niJ n (Ya)H^(ya) (A-l) 

M + (o)M + (-o) = uiJ^yajH^'tYa) (A-2) 


where y - A 2 - c 2 (Ef narsson et al, Equations (5-33), (5-34)). 

The exact function can be evaluated by numerical integration [Mittra 
and Lee, (1971), p. 107, Eqs. (9.8)]. The approximate evaluation 
of these functions proceeds in two different regions, for which the 
Bessel functions are approximated differently. Section A-l deals 
with exact numerical integration. Section A-2 deals with low-frequency 
approximations. Section A-3 deals with high-frequency approximations. 

A-l. Exact Numerical Integration 


The exact expression for L + (a), is given by 

, , ® +ic JlnfuiJ (a/k 2 -z 2 )H^(a/k 2 -z 2 )]dz 

- e*p| JiT / . - 1 " 

' -®MC 


l 


z-a 


(A-3) 


where -Im(k) < c < Im(a) < Im(k). 


Before proceeding further, however, we note that it would be very 
beneficial to not have to integrate out to ®. We note that, for 
Re(z)>Re(k), the arguments of the Bessel functions become essentially 
imaginary. From Abramowitz, Equations (9.6.3) and (9.6.4) 
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J n (ia/z 2 -k 2 ) = e~ 2_ I n (a/z 2 -k 2 ) 


C A-5) 


where K p ( ) and I p ( ) are modified Bessel Functions. 


Taking asymptotic forms (Abramowitz Equations (9.7.1) and (9.7.2)) 
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(A-6) 


We are integrating ln(a/z 2 -k 2 ) from k to °°, but would prefer to inte 
grate ln(l) from k to <», since it is trivial. Hence, we form the 
equation 


-ia 



L + (a)L + (-a) 




(A-7a) 


/-ia(k+a)L + (a) /-ia(k-a)L + (-a) = naJ^Ya^ 1 ^ (ya) (A-7b) 


We can evaluate this by the same integral, except that the argument 
of the logarithm becomes very close to 1 for Re(z)>Re(k). Taking 
c*0 


/-ia(k+a)L + (a)-expL^ / 
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dz 


(A-8) 
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Similarly 


/-ia(k+a) M + (a)=exp 
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> dz 
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(A-9) 


These integrals have the advantage that we can essentially integrate 
on the range where -(ka+C) < z < (ka+C) and the total error thus 
introduced can be made negligible by proper choice of C. 

The complex value of k requires some thought. The argument 
introducing this imaginary part is given by [Einarsson (1966), 
p. 147]. The correct answer is defined as the one obtained in the 
limit where Im(k) -*■ 0. For the purpose of numerical integration, 
two routes present themselves. We can either allow k and a to have 
small but finite imaginary parts, or we can attempt to integrate 
in a symmetric fashion about the singularities at z=a and z=k, when 
they lie on the path of integration, and add in the appropriate 
residues. In fact, the first alternative was chosen. It was evident 
that the factorization functions must be smooth, hence slowly varying 
as functions of complex k and a. Also, experimentation with values 
of loss tangent indicated the results were relatively insensitive 
to the loss tangent. Hence, for purposes of numerical integration 
complex values of k and a (being k and a ) are defined by 


k Q = (1 + i.005)k (A-10) 

a Q = a + i.0025 k (A-11) 

The results of numerical integration for values of n (the order of 
the Bessel functions listed as n in Table A-1) ranging from 0 to 
3 are presented in Table A-l. It was found, by way of confirmation, 
that there was excellent agreement with both DC and asymptotic forms. 
The $M + function is defined by 


$M + ( a) = a(k + a)M + ( a) 


(A-12) 


It remains finite as k -*■ 0. 
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Table A-l 

VALUES OF FACTORIZATION FUNCTIONS L+ AND M+ COMPUTED BY 
NUMERICAL INTEGRATION 
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A-2. Low-Frequency Approximations 

Derivation of 'DC' forms for the factorization functions is 
not difficult. They can be obtained directly from the small argument 
approximations for Bessel and Hankel functions. For n ^ 0, we get 
(Abramowitz, Equations (9.1.7) and (9.1.9)) 


L+(a)L + (-a) 

= niJ^yajH^^Ya) 

(A-13a) 



(A-13b) 


^ n 

(A-l3c) 

L + (a) 

/n 


(A-14) 

M + (a)M + (-a) 

= 7 fiJ^(Y a )H^ 1} ' (ya) 

(A-15a) 


i a n-1 

i ?( ^ } 1 n! r n 

~Tn^ryr ■ 

(A-l5b) 


M + (a) 

For n = 0, we note that 

J’(z) = - ^U) 

H^ }, (z) = - H^(z) 


(A-15c) 

(A—16) 


(A-17a) 

(A-17b) 
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Hence, 

M + (n=0,k,a) = L + (n=l,k,a) (A-18) 

From Abramowitz Equations (9.1.7) and (9.1.8), for n = 0 

L + (a)L + (-ot) = mi J Q (Ya)H^ (ya) (A-19a) 

if ni-1 * — Jtn(ya) £ - 2Zr\ ya (A-19b) 

$ -ln(a 2 (k 2 -a 2 )) (A-19c) 

The factorization of this DC form is non-trivial. Based on[Weinstein 

(1969), p. 335-36] we may approximate L (a) ;f i/2 — f or 

a « k « 1. An(2ka) 

Having thus obtained the first low-frequency approximation 
with ease, it is extraordinarily difficult to obtain a better approxi¬ 
mation analytically. Lee, Jamnejad, and Mittra present a simple 
derivation. Unfortunately it does not agree with the results of 
numerical integration. This can be seen by comparison of their 
formulas with the plot in Figure A-l, which shows the trajectories 
of L + (k) and M + (k) for n=l computed by numerical integration. The 
inclusion of additional terms from the small argument approximation 
does not improve matters. The approach taken, therefore, is that 
for small values of ka, the values of L + and M + are computed by 
numerical integration, and interpolation is used. 

It turns out that, for L + (n=l,a=k), the values over a surprisingly 
large range (0 < ka < 2) can be approximated by 


L 


+ 


1-i.0584 ka 
(l+i.68 ka) 


(A-20) 


However, this does not admit any generalization, nor does it possess 
any theoretical justification. 


High Frequence 


itotic Approximations 


Using the large argument approximations for Bessel and Hankel 
functions, we obtain (from Abramowitz Equations (9.2.17) through 
(9.2.20)) 
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F igure A-l. Factorization functions L+, M+, $M+ for n=1,a=k. 
Arrows indicate increasing k. 
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(A-21a) 


L + (a)L + (-ot) = niJ n (Ya)Hj | 1) (ya) 


i®(Ya) 

=iriM n (Ya)cos9 n (Ya)M n (Ya)e 

(A-21b) 

TTi 2 i2 0 n (Ya) 

* | 1 (M n (Ya)) 2 (l+e n ) 

(A-21c) 

= HJ^YajH^'fYa) 

(A-22a) 

1* n (Ya) 

= TTiN n (Ya)cos<)) n (Ya)N n (Ya)e 

(A-22b) 

Tti 2 12* (Ya) 

= J i (N n (Ya)) 2 (l+e n ) • 

(A-22c) 


Where, from Abramowitz, Equations (9.2.28) through (9.2.31) 


V 2 ) 2 'k 


i + 


4n 2 - 


8z 




0 n (z) 'I - ({[ + 1)* + 


"5T 


Wrf'TS 


4n 2 -3 

8z 2 


V2 ) , 2 . ( n. l )ir+ 4r^3 


(A-23) 

(A-24) 

(A-25) 

(A-26) 


For L + (a) we observe that we can separately factor the two expres¬ 
sions, making L + (a) the product of two factors, both of which are 
analytic in the upper half-plane. 


L + (a) = Lj(a)[_j I (oO 


(A-27) 
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(A-31) 


M {(«) ' 


e iir/4 
J a(k+ot) 


/ 


k 2 a 2 - 4 - n - - J r- + aa 


a(k+a) 


In evaluating L^(ot) and M^(a), we follow the derivation of 
et al, (1975) . 


L^WL^-a) = 1 + e n 


129 (ya) 


M| I ( Qt )M^ I (-a) = 1 + e n 


12<t> (Ya) 


II, 


The exact expression for L* (a) is given by 


L^a) 


= exp< 


_ CP 

1 r £n 
2rf J 


1 + e 


129 n (ya) 


z-a 


- dz 


Let 


,. j»(,, ; 29 " (Ta) ) 


z-a 


dz 


Then the identity 


00 ( n nH ‘ 1 v m 

ln(l + x) = l ; |x| < 1 

m=l 


can be used to obtain 


i2me„(va) 

® , ■, xnH-1 ® „ n v ' 


m=l 


z-a 


dz 


provided Im(e n ) > 0. 


Mittra 

(A-32a) 

(A-32b) 

(A-33) 

(A-34) 

(A-35) 

(A-36) 
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Let 


I 


m 


=; 


12m6 (ya) 

e n dz 


z-a 


(A-37) 


where 


e„(ra) 

2 

* , / n , 1 \ L 4n -1 
' Ya " + T )lT + ~5^a 

(A-38) 

ya = a' 


(A—39) 

z = k 

sinx 

(A-40a) 

dz = k 

cosx dx 

(A-40b) 

ya = ^ 

/ 2 2 2 

vk -k sin x = ka cosx 

(A-40c) 


The contours of integration in the z and t plane are shown in Figures 
A-2 and A-3, respectively. The contour in the x plane is selected 
so that ya has the right sign, and the integrand vanishes at both 
ends of the contour. The contour in the x plane is then deformed 
to the steepest descent contour. Since no poles are passed over 
in this deformation (for Re(o)>0) the integral is unchanged. Hence 


I 


m 


f 

C 


i2m0 n (kacosT) 

ksint-a 


kcosx dx 


- f 

C SDP 


cosx 

sinr-f 


i2m6 (kacosx) 
e n dx 


With 0 n (kacosT) = kacosx - (£ + ^)tt + 


4n 2 -l 

8kacosx 


(A—41) 


(A—42) 


. 4n 2 - 
1 8(ka) 2 


(A-43) 
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I m is an integral in the form 

f F(t) e 1 ^ 1 * dx 
C SDP 


where F(x) = — ° - ST - 
sinx 


a 

k 


k = ka 
f (x) = 


i2m9 n (kacosx) 

fa 


= i 2 m(cosx - + 3 -) 


X L 

+ -L-) 

cosx 


(A-44a) 

(A-44b) 


(A-44c) 


This can be readily evaluated by asymptotic techniques. We make 
a first approximation by ignoring the pole near the saddle point. 


X. sinx 

f' (x ) =0 = i 2 m (-sinx +-*—) 

b COS X 


(A-45) 


= 0 


(A-46) 


Then f"(x s ) - i2m(-cosx + 2 " C 3 S T - X L ) 


COS X 


= -i2m(l-X L ) 


(A-47) 


The asymptotic value of the integral is then 


«c f (T s ) -r 

! m " e 


/iw 

V<^TU 


F( t $ ) 


where ^ is the angle of the C SDp through the saddle point. 


(A-48) 












i2me n (ka) 


=-TcT 


F(t s ) = " l 


f"(T $ ) = -i2m(l-X L ) 


(A-49a) 


(A-49b) 


(A-49c) 


Thence 


I_ e 
m 


i2m0 n (ka) e " iir/4 


(- -) 
a 


(A-50) 


,. r ( .D" e i;me " (ta) g ' W4 . 

m a' 

m= I 


k 

j' — e 
a 


- i tt/4 


Ji 


[-ir e 


i2m0 n (ka) 


(A-51) 


The indicated sum is shown by [Chuang, Liang, and Lee (1975), p. 773] 
to be a modified Lerch function. They present a development of its 
properties and a transformation to simplify its evaluation. 


* (-l) m e ,2m6 n (ka) - e i2™ P 

A. —^T75- = 1 T77 


1 372- = LCp.3/2) 

m= i m 


(A-52) 


where + irmr + i2m0 n (ka) = i2ump 


or p = 


e„(ka) 


(A-53) 


This function has two useful and simple properties. Since the 
parameter p appears only in the exponent, and 
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e i2wnp+i2np _ e i2imip 
L(p+1,3/2) = L(p,3/2) 


(A-54) 


Furthermore e ^2mn(l-p) = e " i ^ 1,rnp 

L(l-p,3/2) = L*(p,3/2), where * denotes the complex conjugate. 

(A-55) 

By using these relations, it is possible to transform any value of 
p to the range 0 < p < 0.5. In this range, the series solution is 
extremely accurate 


L(p,v) 


r(1-',)e i(1 - v>,,/2 , r i m c(v-in)(2np) m 
(2np) ~ v mio 


(A-56) 


For v=3/2, this gives 


L(P,3/2) = -2ir ^-(l-i) + l (R + iI n P)p 2n (A-57) 

n=0 n n 


where values of R n and I are given by Table I of Chuang, Liang, 
and Lee (Note sign errorin Equation (35) of Chuang et al (1975) 
corrected here). 


We can similarly find an asymptotic form for M (a). The analysis 

+ 2 

is identical with 8_ (ka) replaced by <fe (ka) and X, -► X M = -- n 

V(ka) n L M 8(ka) Z 

This results in p = —- + 1/2. Otherwise the expression is 

unchanged. " 


' L»(a) ' 


*l,m 

""(a) 

, s exp < 



l J 


(A-58) 


Where 


I 


L 


k - i it/ 4 
a e 




ka(1-X L ) 


L(p,3/2) 


9 (ka) 

p * -Ajj— + 1/2 (A-59a) 
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>M ' 3 e ' in/4 L(M/2) P ’ ^ - VZ (A - 59b) 

A comparison of these results with those obtained by numerical inte¬ 
gration is instructive. For a=k, the agreement is very good down 
to the first zero of J n (ka) or J^ka) (for L + and M + respectively) 

if the proper form of 0 n (ka) or 4> n (ka) is chosen. This is illustrated 

in Figure A-4 for M+(n=l, a=k); the exact value of M + is compared 
to that obtained by using the two term and three term approximations 
to <Mka). It can be seen that a rather sharp minimum occurs for 
the exact form of M + for that value of ka for which J,'(ka)=0. For 
the two-term and three-term forms of <Mka), the minimum occurs at 
the value for which cos<^(ka)=0. Thus, the more accurate 9 n or 4> n 
should, in general, lead to a more accurate approximation for L 
or M + . 


L + and M + are also functions of a/k; we must explore the ac 
curacy of this approximation when this ratio varies. This is 
tant since a/< = cos9 for scattering, radiation, and coupling 
. 

and aj k 




impor- 

problems. 


) for coupling and radiation of waveguide modes. 


In Figure A-£, the value of L + (n=l) is plotted in the complex plane, 
where the main curve traces out the real and imaginary parts of L + 
with *=k and D/X as a parameter. The dashed curves then indicate 
the trajectories which occur when varying a/k while holding ka (or 
D/X) fixed. It is immediately evident that the preceding approxi¬ 
mation does not, in general, capture the nature of this behavior. 

The reason for this can be seen from Figure A-3. As a/k varies, 
the pole located at T p = sin (a/k) moves closer to the saddle point 
at t s =0. We totally ignored this pole originally. Without belaboring 
the point, it will only be stated that rederiving the approximation 
incorporating the effect of the pole produced only a marginal im¬ 
provement. It was therefore necessary to include two terms in the 
asymptotic expansion. This derivation follows. 


With the same definitions of K,f( t) and F(x), we have the 
asymptotic approximation 
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Figure A-4. Factorization function M+ for n=l, ot=k 
calculated by numerical integration and asymptotic 
approximation based on two and three term 
expressions for phase of 
Bessel function. 
















». Factorization function L+ for n=l. Solid curve gives a=k, varying k 
curves vary a/k holding D/X fixed to the values specified by arrows. 

□ = exact value (by numerical integration) 
x = first asymptotic approximation. 








Kf(Tj i<J> 

I_ * e e 

m 


/ 


2*_ 

<f"(T c 


F(x s ) 


+ 1A (1 - f(KA)) 


i F <V 


3f {1V) (T s )f-(T s )-5f-'(T s ) 2 

3f"(T $ ) 3 


+ F'(t s ) 


f ,M (T s ) 

f"(T s P 



(A-60) 


where f(X) is the transition function defined and discussed by 

[Kouyoumjian and Pathak (1974), p. 1453]. A is a distance parameter, 

relating to the separation of x and x , defined by 

s p 


A = i(f(x $ ) - f(x p )) (A-61) 

The saddle point is unchanged, so x s =0. The pole is located at 
T p - Sin’’ 2 

We thus derive Table A-2. 


(A— 62 ) 












TABLE A-2 

SUMMARY OF ASYMPTOTIC FUNCTIONS AT SADDLE POINT IN t-PLANE 



Function form 

Value at 

f(T) 

i2m(cosx-( ? + + cQST ) 

i2m(l-(£ + 

f'(T) 

i2m(-sinT+X. s1n I-) 

L COS i 

0 

f(T) 

2 

i2m(-cosT+X, 2 ~-|- T ) 

L cos J x 

i2m(-l+X L ) 

f"'(x) 

i2m(sint+X. - sin I (6-cos 2 t)) 

L COS^T 

0 

f (lV) (T) 

i2m(cosT+x. co 5 4 t-20cos 2 t*24, 

L cos 5 t 

i2m(l+5X L ) 

F(T) 

COST 

sinx-a/k 

-k/a 

F' (t) 

a/k sinx-I 
(sim-a/k) 2 

1 

Qj*" 

ro 

F"(t) 

2 

-a/k sinxcosx-(a/k) cosx+2cosx 

i - 2 A 3 

a 'a 

(sinx-a/k) 2 



























„ k / 7T 

H/ka(i-x L ; 


e -in/4 .-1/2 e i2m9 n < ka ) 


1 + 


1 1+ 5X l 1-2$) 


2 - 


L H + \> - 1+ \ J 


= m6 


= 2 


/ 


l-(o/k) 2 - 1 + X L ( 1 


& 


(«/k)‘ 


-3+9X l + 8$) (1-X L ) 
8(1-X l ) 2 ~ 


_ k / if 
- 


-i-rr/4 


A, 12me (ka) f } 

j'-- e y + A Q 6(l-f (mka6)) > 

vfiT ' J 


* / . s m+l 

* i 


m=l 


m m 


vr m= 


nT A 


oo , n n*l (-A,) 12me (ka) f 'I 

I, ^-4- • " ('*A oS (l- ? (m kaS ))J 


® , n m 12m0 (ka) f 'l 

1 l V2 e <l+A 0 «(l-f(mka5)) 

m=l r ^ J 

' Aj |L(p,3/2)+A 0 5G(p,q)j 


(A-65b) 

(A-66) 

(A—67) 

(A-68) 

(A—69) 


(A-70) 
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where 




fi(p.q) = I m- 3/2 V' 2 ™Pn-lf(mq)) (A-71) 

m=1 

6_(ka) 

p = - + 1/2 (A-72a) 

q = kafi (A-72b) 

We could follow exactly the same analysis to obtain M+, except that 
9 n ( ka ) -*■ 4 > n (ka) 


m 


4n 2 +3 

8(ka) 2 


(A-73) 


These formulas produce a very good agreement with the results of 
numerical integration, down to a/k=0.2. 















s 


- 2i 

4>cJ> ” k(1+cose i )(l+cose s ) 


l Q e n cos n 4>s 


J^(kasin0.) 
M + (kcos0 1 -) 


J^(kasin0 s ) 

M + (kcos0 s ) 


(l+cos0 i )(l+cose s ) fjj 

2(cos@ i +cos0 s ) + 7^2 
1 11 1 n - 


(B-4) 


, . 41 ? . J n (kasin9 s) j ;< kastrl9 l ) f n 

^ 04 ) k(l+cos0 i ) sin n<p s sin0 s L + (kcos0 s ) M + (kcos6 i ) * 

(B-5) 


_ 4i ■ ^(kas1ne s ) J^kasine,) f„ 

b 4>0 ” k(l+cos0 s ) n fi sin n *s M + (kcos0 s ) sin0 i L + (kcos0 i ) ^2 

1 n 

(B-6) 

Where 

4» i = it is assumed 
nL + (k) 

f n = 2k5H + 0c) (B_7) 


1 


e 


n 



n = 0 

n = 1,2,3... 


(B-8) 


It was found that Chuang et al (1975) contained several misprints. 
A correction letter is reproduced on the following page. 




111 




AP 11 -44 


Cal. I 


Coneclion to “High Frequency Scattering front an Open- 
Ended Semi-Infinite Cylinder” 

C. A.CHUANG, MEMBER. IEEE,CHARLES S. LIANG, 
MEMBER, IEEE, AND SHUNG-WU LEE, 

SENIOR MEMBER, IEEE 

In the above paper', there were six sign misprints. 

1) Time factor should have read exp( -tu>r). 

2) In (1), x should have read (-jr). 

3) In (4), (1 + cos 6) in numeralor should have read (1 - 
cos 8). 

4) In (7), (cos 8 + cos 8 0 ) should have read (-l)(cos 6 + 
cos 8 0 ). ■ 

5) In (11), (a - o’) should have read (a' - a). 

6) In (35), (I + i) should have read (1 - r). 

the above sign misprints do not affect any other 
equations or numerical results. 

The authors wish to thank T. W. Johnson and D. L. Moffatt 
for bringing some of these errors to their attention. 

C. A. Chuang is with Aeronutronic-Ford/WDL, Palo Alto, CA. He 
is now at 552 Solitaire Dr., Indialantic, FL 32903. 

C. S. Liang is w ith the Fort Worth Division, General Dynamics, Fort 
Worth, TX 76101. 

S. W. Lee is with the Department of Electrical Engineering, Univer¬ 
sity of Illinois, Urbana, 1L 61801. 

* C. A. Chuang, C. S. Liang, and S. W. Lee, IEEE Trans. Antenna 
Propagat., vol. AP-23, pp. 770-776, Nov. 1975. 
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B-2. Coupling of Incident Field to Waveguide Modes 


With the same definition of the incident field (as in Equation 
(B-2)) the axial fields within the waveguide can be represented as 
follows: 


E 


z 



00 


l 

m=l 


K" cos n * | e - 1c w 

A^Lsin r>4> i ^n^nrn^ 

[ m J 


(B-9) 



oo oo 
1 

n=0 m= 


y 0 

B„ m sin n$ 
nm 

B^ cos n<f> 
nm 




•iot' z 
nm 


(B-10) 


^n^ J 'nm 


)=0 


J:(j' )=0 
n w nm _ 


’nm 


nm 



The top line of Equation (B-9) indicates that Eq produces 
E z with cos n<p dependence 

The bottom line of Equation (B-9) indicates that produces 
E z with sin n<j> dependence v 


e L,(°t_m) J n (kasin9.) 

l - n nm + nm n v 

nm a a L,(kcosQ,) kasine. 

nm + i i 


f n + <°^)n-c°s6. ! 

7Tj; 2 <‘W kcos6 i> 


(B-ll) 
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(B—12) 


A* , . 3 m 
nm 


- h+{> L. j;o<«ine,) ,JL > 

ka a 2 (l+cos0.) M + (kcos0 i ) n 1 l-f‘ 

nm ' 1 ' r 


e , 

B nm = *- J~ 

■W 1 ' 77T> 

J nm 


M.(a' ) J_(kasin@,) f n 
+ _ nm n' i ' n 

L + (kcos0 i ) kasi‘n0 i ~^c 


(B-13) 


J* = 
nm 


2e i n (k+a' ) 
n ' nm 7 


Km** 1 ' T77 )(1+cos0 i ) 

J nm 


M + («* ) 

n jkgrse,) Ontkasine,) . 


f„ (k.a^)(Hcos0.) 
Tf! 2(kcose,-a^) 


(B-14) 


B-3. Radiation Patterns From an Open Waveguide 

The incident field is assumed to be a single waveguide mode, 
either TE or TM, with axial fields given by 


r J n»™ 0/a >_ 1( V 

2 ‘ ™ TOJ- cos "* e 

n' rim' 
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The radiated field is given by 
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B-4. Internal Reflection and Mode Conversion 


field 


With incident field as described in Section B-3, the scattered 
is given by 
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APPENDIX C 

ELEMENTS OF DYADIC GREEN'S FUNCTION 
FOR A CIRCULAR WAVEGUIDE 


The dyadic electric Green's function, omitting the source- 
singular term, is given by Tai (1971). It is assumed, in this appen¬ 
dix, that the azimuthal dependence for Ep, E , Op, J , and H* is 
cos n$, and for E*, 0*, H p , and H is sin n<ju This azimuthal depen¬ 
dence is understood and suppressed. The summation over n is also 
suppressed. 
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